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ABSTRACT 

We present a new analysis of a Keck spectrum of Q00 14+813, with a new 
model of the z = 3.32 absorbing cloud. We fit Lyman series absorption of the 
dominant HI components, including components identified with hydrogen and 
associated deuterium at two discrete velocities, with thermal line broadening of 
each species. The deuterium features are too narrow to be interlopers, and the 
good agreement in temperature and redshift with their hydrogen counterparts 
confirms the identification as deuterium. The abundance is measured to be D/H 
= 1.9 ± 0.5 and 1.9 ± 0.4 x 10 -4 in the two components, with an independent 
lower limit of D/H > 1.3 x 10~ 4 for the sum, derived from the Lyman Limit 
opacity. The impact on cosmological theory is briefly discussed. 

Subject headings: cosmology: observations — quasars: individual (0014+813) 
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1. Confirming Deuterium 

The hydrogen cloud at redshift z = 3.32 along the line of sight to the quasar 
Q00 14+813, shown by Chaffee et al (1985,86) to be uniquely well suited for measuring 
primordial abundances, was found recently to show an absorption feature at the precise 
wavelength predicted for cosmic deuterium flSongaila, Cowie, Hogau and Rugers 19941 , 



hereafter SCHR, and Carswell et al. 1994| ). The implied D/H abundance ratio from this 
feature, although it fits very well into the Big Bang predictions for primordial abundances 
( |CJopi et al. 1995a| , par 1995| , [Fields and Ulive 1996] , |Ulive et al. 1996] ), is much higher than 



that usually derived from models of Galactic chemical evolution normalized by solar system 
3 He abundances ( |Galli et al. 1995| , [Wilson and Rood 1994| , |Hata et al. 1995| , popi et aT. 



1995b|) , and also much higher than the best estimate of D/H in another QSO absorber, on 



the line of sight to Q1937-1009 QTytler and Fan 1994|) . The simplest explanation of this 



discrepancy is that the feature in Q0014+0813 is not caused by deuterium at all, but is 
simply a low column density hydrogen cloud that happens to lie at a somewhat different 
velocity from the bulk of the local hydrogen, at the redshift expected for deuterium (SCHR, 
and |Stcigman 1994|). Although the redshift coincidence is quite precise (±5 km/sec), the 



chance of such a coincidence is still not negligible, and was estimated to be at least a few 
percent, ignoring possible correlations. 

We can eliminate the possibility of interlopers even in a single system if the candidate 
D lines are very narrow. The linewidth, characterized by the Doppler parameter in a 
thermal profile fit, b = \J^- is 13T 4 1//2 km s -1 for hydrogen, 9T 4 1//2 km s -1 for deuterium. 
The potential interloper population is observed to have lines below b = 20km s _1 only very 
rarely (|Hu et al 1995|) , so a narrower line is much more naturally identified with deuterium 
than with a chance interloper. Also in very narrow lines, thermal broadening is likely to 
dominate the linewidth, so the ratio of the fitted Doppler parameter for D candidates to 
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their H counterparts is an important clue. Here we use evidence from line shapes to confirm 
the identification of the deuterium features and sharpen the estimate of the primordial 
deuterium abundance. 

2. Observations and Reduction 

We use the same data described in SCHR: six exposures of 40 minutes each were 
obtained of Q0014+813, at the Keck telescope on November 11, 1993, using the HIRES 
echelle spectrograph with a 1.2" x 3.8" slit. A resolution R=36,000 was obtained over the 
covered wavelength range of 3500 to 6080 A. Unlike SCHR, the raw data were reduced using 
IRAF. The CCD frames were overscan corrected, bias subtracted and flat fielded using 
CCDPROC. After verifying the registration of the frames, the frames were then stacked 
using IMSUM, with rejection of the highest pixel value among the frames, in order to reject 
cosmic rays. Fewer cosmic rays were left using this method than by creating a median with 
IMCOM. The stacked images were then used for extraction of the individual orders of the 
spectrum using DOECSLIT, where the standard star frame was used to trace the orders on 
the CCD and ThAr calibration lamp data were used to obtain the wavelength calibration. 
No attempts were made to flux calibrate the spectrum. The resulting extracted, dispersion 
corrected echelle spectrum was then run through RVCORRECT to obtain the correction 
for radial velocity of the observer with respect to the QSO, so redshifts quoted here are in 
vacuo, relative to the local standard of rest. The data were not smoothed, as they were 
in SCHR; the new results presented here are due primarily to using the full instrumental 
resolution of the raw data. 
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3. Model Fits for HI and DI 

After examining the Lyman series of this absorber complex, we decided to use 
Lyman-a, (3, 7, 5, (, and k through 17 for establishing fits. Lyman-e was thrown out due 
to the presence of significant absorption from Lya- forest interlopers. A value for log Nhi of 
17.3 (estimated from the Lyman limit optical depth) was used as an initial guideline for the 
column density of the bluest part of the complex, where most of the HI is located, as is 
shown by the higher order lines. The code used to perform the model fits to the spectrum 
is VPFIT, which models the formation of an absorption spectrum as a series of thermally 
broadened components at discrete velocities ( Carswell et al. 1987] , Webb 1987). For each 



component VPFIT determines the redshift (z), Doppler parameter (b) and column density 
(N) by fitting the data with Voigt profiles convolved with a Gaussian instrument profile, 
and makes formal error estimates from the covariance matrix parameters, for each of the 
calculated z, b and N, based on the reduced x 2 . (Note that although these are "true" la 
errors, including the total fitting uncertainty in each parameter, the probability distribution 
is highly non-Gaussian, so they do not translate directly into "confidence intervals" . The 
results quoted in SCHR were used as a starting point for the fits performed by VPFIT. 

The Doppler parameter for the D-component fitted to the smoothed data by SCHR, 
14 km/s (assuming thermal broadening only), is too large for a good fit to the unsmoothed 
data. A narrower fit is required both by the steep profile edge of the D feature at the blue 
side of this absorber complex, and by a sharp spike in the center, which may at first sight 
appear to be noise. However, since in the top part of this spike there are 3 channels, each 
of which has between 6 and 7 a counts above the bottom of the absorption dip, this is 
statistically very unlikely to be a noise feature. The feature also appears if the same data 
are reduced using the independent software used by SCHR, indicating it is not due to the 
different procedure performed here (A. Songaila, private communication). A much better 
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fit was obtained using two HI absorbers instead of one, each of which has a log N HI of 
around 16.8, both with corresponding DI absorption lines. This two-component fit is also 
preferred by the HI lines on their own, without the split Da feature; by adding the three 
new parameters corresponding to the one new HI component, the reduced x 2 P er degree of 
freedom for the best fit is equal to 0.98, instead of the 6.1 obtained for the best fit with 
only one main HI hydrogen absorber. This is mainly due to the improved fit in the higher 
Lyman series lines, [more here] 

The best overall fit was obtained by the components as presented in Table 1, where the 
components of interest for the deuterium abundance are numbers 11 to 14. Figure 1 shows 
plots of the data and the fit to Ly-a through 5, plus (, all plotted on the same velocity 
scale. Most interloper features in the higher orders were omitted from the fit. Some areas 
of poorly fitted excess flux in 8 and ( must be due to artifacts or noise, as corresponding 
flux does not appear in higher orders. The highest order line fits are shown in Figure 2. 
Although many parameters of the minor components are not well determined by the fit, 
the HI columns of the dominant components (12 and 14) are well constrained by the line 
fits in the unsaturated high Lyman series, and the DI columns are well-determined by 
the unsaturated Lj-a lines, so the absolute deuterium abundances can be estimated fairly 
accurately. 

The D/H ratios for both of the main absorbers are in good agreement with each 
other, as well as with the value obtained by SCHR. For the z = 3.320482 absorber we find 
(D/H) = io- 3 - 73±0 - 12 , and for the z = 3.320790 absorber we find (D/H) = i(r 3 - 72±0 - 09 , 
where we have added the errors in \ogN HI and \ogN DI in quadrature. We therefore 
expect the D/H ratio for these absorbers to lie within the ranges of OtoJ x 10~ 4 and 
1.9lo4 x 10~ 4 , respectively. The total D column is 2.6 ± 0.4 x 10 13 , which also gives 
(D/H) = 1.9 ± 0.4 x 10~ 4 when compared to total hydrogen column in the two main 
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absorbers. These are consistent with the estimates made by SCHR and Carswell et al. 
(1994), using lower resolution. 

The total N HI in the two main absorbers from this fit is 1.37 ±0.16 x 10 17 (the errors 
here obtained by adding in quadrature the errors for each), which agrees with the total 
Nhi of the whole absorber complex estimated from the Lyman break. This can be seen 
qualitatively in Figure 2, from the good match of the fitted absorber model at the Lyman 
Limit to the regions of highest flux shortward of the limit. The flux beyond the Lyman limit 
can be used to derive an upper limit for the total hydrogen column in the absorber complex, 
which yields a fairly stringent lower limit on the D/H ratio independent of HI line fits. In 
order to try and circumvent any leftover instrument effects in this part of the spectrum, the 
continuum of the QSO and the standard star were assumed to have the same shape over 
this narrow range (50 A), and the shape of the standard star continuum was fit to match 
the parts of the QSO spectrum at the red end of the order, in regions without conspicuous 
absorption. This gave an estimate for the continuum level, which was used both for the high 
order line fit and the continuum level beyond the Lyman Limit, and appears at the blue 
end of the fit in figure 2. The level of the signal beyond the Lyman Limit was estimated in 
two ways. The simplest of these used the mean of the counts per channel, blueward of the 
Lyman Limit in this order. This gives an upper limit to the opacity of the cloud complex, 
and hence an upper limit to the total hydrogen column of the cloud, which can be used to 
find a very conservative lower limit to the D/H ratio of the two absorbers. The second uses 
the mean of only those regions blueward of the Lyman Limit where the counts per channel 
are at least 1 a above the zero-level, over a range of at least 1 A. Two examples are the 
regions around 3932 A and 3941 A, each of which is a few A wide. Since additional Ly-a 
and Ly-limit absorption must still occur in these areas, the value for the opacity determined 
from this estimate of the Lyman Limit signal is also an upper limit. The two methods yield 
upper limits to the total HI column density of 2.4 x 10 17 and 1.7 x 10 17 cm~ 2 , respectively. 
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This limit almost exactly matches the sum of all the fitted components, (including 
in particular number 16 in table 1 in addition to the two main components), so there 
is almost no room to hide more than a small fraction of additional HI anywhere in this 
complex beyond that already accounted for in the fits. Taking the lower limit on D to be 
2.2 x 10 13 cm~ 2 , which is well constrained by the unsaturated Da, therefore yields a firm 
lower limit (D/H) > 1.3 x 10 -4 , in good accordance with the lower limits determined from 
the hydrogen line fits. Therefore, if the D candidate lines are indeed caused by deuterium, 
D/H must be high. 

Most significantly for confirmation of deuterium, the DI lines are narrow, making the 
interloper possibility very unlikely. The fits for the two components are 7.5 ± 1.2 km/s and 
8.8 ±1.1 km/s. Hu et al (1995) find that in 3 QSOs (including this one, using the same 
data we are using, and others with similar resolution and signal-to- noise), only 11 out of 
670 lines are found with b < 14. This corresponds to less than 2%, or only a few lines 
per QSO, in the 6000 resolution elements between Lya and Ly/3, almost all of which are 
identified as metal lines. Since interlopers must be metal lines, any coincidence must be 
purely accidecntal, (not affected by physical Lya autocorrelation), so we can compute the 
interloper probability from Poisson statistics. The probability of a chance coincidence of 
one of these lines at the precise velocity expected for deuterium is then much less than than 
10~ 3 for each component. 

Even more suggestive, a thermal model, likely to be appropriate for such narrow lines, 
predicts that the HI counterparts of the deuterium lines should have b = 10.2 ± 1.7km s _1 
and 12.4 ± 1.6km s _1 , in agreement with the observed values of b = 10.1 ± 2.0km s -1 and 
12.7 ± 0.8km s -1 . This coincidence also argues against random interlopers. 

In fact the HI lines are unusually narrow for forest lines, but this is not surprising 
since this system was selected for high column density, which often goes with low ionization 
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parameter and temperature. We also examined the data to look for metal absorption lines, 
but found only upper limits for CII, CIII, Sill and Silll, all species with several lines in 
regions of the spectrum which are unobstructed by interlopers, and which have reasonable 
S/N. These results are nevertheless consistent with a sensible physical model of the main 
absorbing clouds. Using the models of Donahue and Shull (1991), we adopt a conservative 
limit on the ionization parameter for the gas of log?7 ~ —4 (where U = n^/nn), which 
gives an equilibrium temperature limit T 4 < 1.2 consistent with the narrow width of the H 
and D features. This gives a neutral fraction of the hydrogen of 10 2 , and estimated 

C and Si abundances ~10~ 3 solar. With such a small ionization parameter (U « 10~ 2 ), 
charge-exchange reactions D + + H # D + H + guarantee that H and D are locked to the 
same fractional ionization (J. Black, private communication), so N DI /N HI indeed gives the 
true deuterium abundance of the gas. 

4. Discussion 

These results reinforce the interpretation of the deuterium features in Q0014 as 
deuterium associated with the hydrogen. First, instead of one component we now have two 
cases where the H and D velocities agree. A second coincidence is that both components 
give the same abundance - that is, the estimated DI columns of the individual components 
are in the same ratio as the columns of HI. The third coincidence is that the width of both 
DI features agrees with that predicted for DI at the same temperature as the hydrogen. 
Finally and most significantly, if the features are interlopers they are too narrow to be part 
of the numerous Lja forest, and therefore need to be two of the much rarer narrow and 
uncorrelated metal line population, making the probability of two chance interlopers at 
the precise required velocities vanishingly small. By contrast, the deuterium interpretation 
makes good physical sense, and provides a natural explanation of all of these coinicidences. 
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The low metal abundance makes it seem unlikely that more than a small fraction of the 
gas has been processed through stars, so that a negligible fraction of the initial deuterium 
has been destroyed. Since the big bang is the only known source of deuterium ( |Reeves| 
et al. 19731) , we may take the absorber abundance as our best estimate of the primordial 



abundance, (D/H) p = 1.9 ± 0.4 x 10~ 4 . Standard Big Bang Nucleosynthesis then gives a 
baryon to photon ratio rj = 1.7 ± 0.2 x 10~ 10 , corresponding to Q^h 2 = 6.2 ± 0.8 x 10~ 3 . 
For this value of 77, SBBN predictions are consistent with estimates of cosmic 4 He and 
Li abundances (Copi et al. 1995ab, [Hata et al. 1995 , Par 1995 , Fields and Olive 1996 ). 



The total baryon density is then only about three times larger than the cosmic density of 
known baryons in gas and stellar populations (Persic and Salucci 1992|) , allowing a tidy 
picture where most cosmic baryons reside in or near galaxies and clusters of galaxies rather 
than in a dominant, diffuse intergalactic medium ([Fukugita et al. 1995|) . There are not 



enough extra baryons to make most of the dark matter in galaxy halos, consistent with the 
constraints from MACHO observations ( Alcock et al. 1995 ), and implying that the bulk of 



galactic dark matter is nonbaryonic. 

Our estimated primordial deuterium abundance is much larger than that of the Milky 
Way ( Linsky et al. 1995| ), probably due to burning in stars, which converts it to 3 He. 



The low value of 3 He inferred in the presolar nebula and found in parts of the interstellar 
medium ( |Wilson and Rood 1994] ) can be explained as further conversion of the 3 He to 
heavier elements in stars flHogan 1995| , |Wasserburg et al. 1995| , |Charbonnel 1995] , |Weiss et 



al 1996j , piive et al. 1996]) , rather than as a low primordial value of deuterium. 
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Table 1. VPFIT results for the absorbers in Figure 1. 
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Figures 

Figure 1. Spectral fits to Lyman a, j3, 7, 5, and (, plotted on the same velocity scale. 
Ticks mark the centroids of individual components, with parameters tabulated in Table 
1; the overall model fit comprising the sum of this absorption is the solid line. The ticks 
corresponding to the bulk of the hydrogen column, and their deuterium counterparts, are 
systems numbered 12/14 and 11/13, respectively. [Postscript and GIF figures are available 
at |http://www. astro. Washington. edu/rugers/res.html|| 



Figure 2. Spectral fit of the same absorption model to Lyman 10 through 17. The data 
in the plot (but not the fit) has been boxcar averaged at 7 km/sec to reduce noise. The 
plotted model does not include photoelectric continuum absorption and recovers beyond 
the limit to the QSO continuum; the predicted continuum from the line absorption model 
is the asymptotic value of the line fit curve longwards of the break, which agrees with the 
observed flux beyond the limit. 



